Future CMB Constraints on Early, Cold, or Stressed Dark Energy 
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We investigate future constraints on early dark energy (EDE) achievable by the Planck and 
CMBPol experiments, including cosmic microwave background (CMB) lensing. For the dark energy, 
we include the possibility of clustering through a sound speed cj? < 1 (cold dark energy) and 
anisotropic stresses parameterized with a viscosity parameter c^ is . We discuss the degeneracies 
between cosmological parameters and EDE parameters. In particular we show that the presence of 
anisotropic stresses in EDE models can substantially undermine the determination of the EDE sound 
speed parameter c%. The constraints on EDE primordial energy density are however unaffected. We 
also calculate the future CMB constraints on neutrino masses and find that they are weakened by 
a factor of 2 when allowing for the presence of EDE, and highly biased if it is incorrectly ignored. 



I. INTRODUCTION 

For about a decade cosmological data from cosmic mi- 
crowave background (CMB) anisotropy experiments ([l[, 
, ) , in combination with complementary results from 
galaxy surveys [1, Q and Type la supernovae @, 0] , sug- 
gest in an unequivocable way that the present energy 
budget of the universe is dominated by an exotic form of 
energy coined dark energy. 

The presence of a cosmological constant term A in Ein- 
stein's equation of General Relativity (GR) is the sim- 
plest explanation for dark energy The Lambda cold dark 
matter scenario (ACDM) is a simple model that consis- 
tently accounts for all observations, and has therefore 
emerged as the standard model of cosmology. Despite 
the simplicity of this concordance model, however, the 
presence of a tiny but nonzero cosmological constant is 
vexing, and is not understood from the point of view of 
fundamental theory (see e.g. Q and references therein). 
Dark energy could therefore be different from a cosmo- 
logical constant, and indeed many diverse models are also 
consistent with the data [914TT1]. 

Within the framework of a non-interacting, minimally 
coupled additional component to the energy density, a 
general dark energy fluid and the cosmological constant 
may differ in two main aspects: the latter behaves as a 
homogeneous fluid with a constant energy density, while 
the former is a non-homogeneous fluid with a time de- 
pendent energy density and pressure. A simple way of 
describing these models is by specifying the equation of 
state w = p/ p, where p and p are the dark energy pres- 
sure and density. The cosmological constant corresponds 
to w = — 1, while a general dark energy fluid may have a 
time dependent equation of state w(a) which is as func- 
tion of the scale factor a(t), so that w ^ — 1 in general. 

Density perturbations in the dark energy component 
could also leave an imprint in cosmological observables, 
while A is purely homogeneous. The clustering proper- 



ties of different dark energy models are usually parame- 
terized by an effective sound speed, defined as the ratio 
between the pressure to density perturbations in the rest 
frame of dark energy; c 2 = Sp/Sp (see, e.g., (l2l - [l3 |). 
Moreover, anisotropic stress can also affect the density 
perturbations. For example, in the case of a relativistic 
component, anisotropic stresses act as a form of viscos- 
ity in the fluid and damp density pertubations. If dark 
energy behaves like a relativistic fluid in the past, then 
the effects of viscosity should also be considered. 

To parameterize viscosity in a dark component one can 
introduce the viscous sound speed c 2 is , which controls 
the relationship between velocity/metric shear and the 
anisotropic stress [H, A value of c 2 is = 1/3, for 

example, is what one expects for a relativistic component, 
where anisotropic stress is present and approximates the 
radiative viscosity of a relativistic fluid. The standard as- 
sumption is that c 2 is = 0, which however cuts the Boltz- 
mann hierarchy of perturbations at the quadrupole, forc- 
ing a perfect fluid solution with only density, velocity and 
(isotropic) pressure perturbations. 

Any indication for perturbations in the dark energy 
fluid would falsify a scenario based on the cosmological 
constant. However, since perturbations become observa- 
tionally unimportant as the equation of state approaches 
the cosmological constant value, w = — 1, to detect them 
one needs some period in cosmic history when w differs 
substantially from —1. Such a deviation in if is con- 
strained at late times by the observations, so we are led 
to consider this at early times, along with a non-negligible 
early dark energy density. 

Such early dark energy can arise in some cases of the 
tracking class of dark energy models (see, e.g., In 
particular, in tracing models the dark energy density is a 
constant fraction of the dominant component, radiation 
or matter. If this fraction is non-negligible, dark energy 
could therefore be appreciable not only in the late uni- 
verse but also at early times. Several models of "early" 
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dark energy (EDE, hereafter) have been proposed (e.g. 
[H,[li| and references therein). 

Our paper is organized as follows: in Section II we dis- 
cuss the EDE model and the behaviour of perturbations. 
In Section III we explain the types of CMB data used 
and the forecast method, including the weak lensing sig- 
nal. In Section IV we present our results, and finally in 
Section V we discuss our conclusions. 



II. EARLY DARK ENERGY 
A. Model 

In [l8[ a parametrization for the dark energy density 
parameter n de (a) and equation of state w(a) has been 
proposed to recognize the important feature of early dark 
energy. In this model f2° and are the current dark 
energy and matter density, respectively, and a flat Uni- 
verse is assumed so fi^+fi^ = 1. The model is described 
by: 
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where f2 is the early dark energy component density, 
constant at high rcdshift, a eq is the scale factor at matter- 
radiation equality, and wq = w(a = 1). In Figure [1] we 
plot fide(a) and w(a), for wq = —1, O c = 0.03 and f2jj c = 
0.7. Note the energy density ^de(a) goes to a nonnegligi- 
ble constant in the past (whereas £Ia(<i = 10~ 3 ) w 10 -9 ). 
The dark energy equation of state w(a) clearly shows 3 
different behaviours: w ~ 1/3 during the radiation dom- 
inated era, w ~ during matter domination and, finally, 
w ~ wq in recent epochs. 

Moreover such an EDE model with constant sound 
speed can behave like barotropic dark energy models (see 
e.g. [1^, and Fig. 2 of [21|). These models have an ex- 
plicit relation determining the pressure as a function of 
energy density that bring advantages to overcome the co- 
incidence problem and to predict a value of wo « — 1 at 
late times, considering purely physical properties rather 
than being adopted as phenomenology. 

Recent analyses have placed constraints on EDE us- 
ing the available cosmological datasets and forecasting 
the discriminatory power of future CMB probes such as 
Planck (see e.g. [2ll424l |). As recently shown in particular, 
the effects of EDE could be important when combining 
CMB data with baryonic acoustic oscillation data [25| . In 
this paper we follow the lines of these recent papers and 
we present a forecast for EDE parameters from the near 
future Planck [H and far future CMBPol [27| experi- 
ments. Our work will improve similar recent analyses in 
several aspects. First, we consider the possibility of per- 
turbations in EDE including an anisotropic stress term 
in EDE, parametrized by a viscosity sound speed c v i s 
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FIG. 1: Behaviour of the early dark energy model in energy 
density (solid black line) and equation of state (dotted blue 
line) as a function of the scale factor. 



(see [l2|). If EDE is following an equation of state of a 
rclativistic fluid, anisotropic stresses can be present and 
change in a substantial way the theoretical predictions 
on the CMB angular spectrum. Secondly, we include 
the CMB weak lensing signal, discussing its importance 
in constraining EDE parameters. Finally we also tested 
our results performing a full Monte Carlo Markov Chain 
(MCMC) on Planck synthetic dataset. 



B. Perturbation theory 

Here we briefly review the perturbations in EDE and 
show theoretical predictions for the CMB anisotropy an- 
gular spectra and for the weak lensing CMB signal. 

In the synchronous gauge, the energy-momentum con- 
servation in the Fourier space gives the following equa- 
tions for the evolution of the density and velocity pertur- 
bations (see [H, (Hi) : 
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where S and 9 are the dark energy density perturbation 
and velocity perturbation, h is the metric perturbation 
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source, and — h/2 — 3i] is the scalar potential of the ten- 
sorial metric perturbations. 

The above equations describe various models of dark 
energy; note that even if w(a) is the same for two models, 
they can differ in the perturbations. For a chosen model 
one can implement these relations in a modified version of 
CAMB [29[ and solve the Einstein-Boltzmann equations. 

III. EFFECTS ON THE CMB 
A. CMB Angular Spectra 

As already discussed in the literature (see e.g. [l3[ and 
fl4^. perturbations in a dark energy component with a 
constant equation of state and a negligible energy com- 
ponent in the early universe (i.e. O e = and w(a) = wq) 
affect the CMB anisotropy only on very large angular 
scales, where cosmic variance dominates. The reason is 
that since in this scenario dark energy contributes appre- 
ciable energy density only at late times and is minimally 
coupled with other energy components, changes in the 
CMB spectra can be only induced by the late Integrated 
Sachs- Wolfe (ISW) component. 

As an example, we plot in Figure [2~] the CMB angular 
spectra for different values of Cg and c^ is : the variation 
is only present on large scales (low multipoles). As al- 
ready discussed in the literature, the feasibility of accu- 
rately measuring one of these parameters is strongly un- 
dermined by the presence of cosmic variance. Moreover, 
the effects of the two parameters are not uncorrelated 
with each other, as we show in Figure [3] Fixing c 2 is = I 
or Cg = 1 makes the angular spectra independent of any 
variation of the other parameter (c 2 or c^ is , respectively). 
If one assumes either c 2 is = 1 (shown in the top panel) , or 

= I (bottom panel), one is maximally suppressing the 
perturbations, giving essentially identical power spectra 
for different values of Cg or c 2 is , respectively. This dis- 
cussion is fully compatible with the results presented in 

The net effect of increasing or c 2 is is higher ISW 
power. This reflects the increased potential decay due 
to dark energy; while dark energy perturbations would 
help preserve the potential, increasing or c^ is reduces 
the dark energy perturbation contribution and so eases 
the decay of the potential. For example, A leads to a 
high ISW power today. The effect can be explained more 
mathematically as follows. The metric perturbation, h, 
is a source term in the density equation ([3]) and tends to 
draw dark energy into overdensities of cold dark matter. 

However, for positive c 2 and/or positive c 2 is , the term 
proportional to 9 dominates (on small enough scales) and 
suppresses perturbations. In the case of positive c 2 , this 
can be seen directly from equation ^ , where the term 
proportional to c 2 implies that the sign of 9 is the same 
as that of 5 so that the contribution to 5 has the oppo- 
site sign of 8, leading to suppression (a comparison of the 
magnitudes of the different terms shows that the suppres- 



sion becomes dominant roughly for scales k > c~ (a /a)). 
Thus S gets smaller when dark energy begins to dominate 
and the ISW effect is enhanced when one increases the 
sound speed. In the case of positive c 2 is , it follows from 
the sign of the metric terms in equation ([5]) that a ends 
up with the same sign as 5, again giving a contribution 
to the 9 equation of the same sign as 8. Therefore, as 
the dark energy becomes dominant, the overall density 
structure is also smaller when c 2 is is larger, and the ISW 
effect is amplified again. 
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multipole 

FIG. 2: Effect of the sound speed (Top Panel) and viscosity 
(Bottom Panel) on the CMB spectrum for Q e — and a 
constant equation of state w — —0.8. 

It is interesting to investigate if this competition be- 
tween Cg and c 2 is is still present in the case of a EDE 
scenario. For dark energy present at early epochs it may 
also contribute to the early Integrated Sachs Wolfe ef- 
fect. In Figure H] we plot the same spectra as in Figure [3] 
but now with an EDE contribution with f2 e = 0.03. We 
see that now the spectra show a small difference around 
the first peak due to the different early integrated Sachs 
Wolfe effect. While the differences are small it is impor- 
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FIG. 3: As cji s (Top Panel set) or (Bottom Panel set) 
approaches 1, the ISW component of the CMB spectrum sat- 
urates, bringing essentially identical power spectra for differ- 
ent values of the other parameter, i.e. or c^ is respectively. 
The bottom half of each set shows the fractional deviation in 
power among models. 



tant to notice that at these scales the cosmic variance is 
significantly smaller than at large scales where the late- 
time ISW effect is important. 

In addition to the (early and late) ISW effect, the pres- 
ence of EDE also affects the evolution of the acoustic 
oscillations before recombination, leading to a signature 
at larger Z's than the ISW. If the sound speeds are in- 



FIG. 4: Same CMB spectra as in Figure [3] but now with an 
EDE component with early energy density Q c = 0.03. 



creased, EDE perturbations get more suppressed, lead- 
ing to a stronger decay of the metric perturbations. This 
in turn leads to a stronger boost of the amplitude of the 
acoustic oscillations. The (subtle) damping in the sec- 
ond peak is a sign that the potentials have not decayed 
as much as when perturbations are unimportant. 

The ISW behaviour is better shown in Figure [5] where 
we plot just the ISW component of the temperature CMB 
anisotropy angular power spectrum. As we can clearly 
see, the behaviour of the ISW angular spectrum can be 
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evidently divided into a contribution from the late ISW 
effect on large angular scales {I < 30) and a contribution 
from the early ISW, producing a peak on degree scales at 
£ ~ 120. While variations on large scales are negligible 
compared to cosmic variance errors, perturbations intro- 
duce signal via the early ISW term that is more signifi- 
cant. We can therefore expect that in the EDE scenario 
perturbations can play a more significant role than in a 
standard late dark energy scenario. The perturbations 
also influence gravitational lensing of the CMB, as we 
discuss in the next section. 



parameters (see e.g. [3(1 HU), since it is strongly con- 
nected with the growth of perturbations and gravita- 
tional potentials. The effect of weak lensing is to remap 
the direction of observation (see e.g. [32h34| ) from n to 
n' = n + d(n) where d(n) is the lensing deflection angle. 

The lensing deflection angle power spectrum, or equiv- 
alently the convergence power spectrum, is related to the 
lensing potential spectrum Cf^, through : 

Cf d = 1(1 + l)Cf+ . (6) 
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FIG. 5: As Figure [4] with an early dark energy density fi e = 
0.03, but now focusing on only the ISW component to show 
the effects of the sound speeds. The left rise is due to the late 
ISW effect while the bump is principally coming from early 
ISW. 



FIG. 6: Lensing potential power spectra for different standard 
and early dark energy scenarios with and cj ia varying from 
to 1. 

Figure [5] shows the lensing potential angular spectra 
for scenarios with and without EDE and for different 
values of and c^ is . The plot shows a nontrivial depen- 
dence of the lensing angular spectrum on c^, c^ is , and 
f2 e , with some degeneracies clearly present. Basically, 
suppressing perturbations by taking c^ is = 1 or = 1 
(or both) are nearly equivalent. Only when perturba- 
tions are maximally allowed, through c^ is = and = 
together, is the lensing power significantly enhanced. In 
this case, early dark energy plays a major role, yielding 
a 30% enhancement in power, while the model with no 
early dark energy only sees a ~ 6% boost relative to its 
no-perturbation case. We therefore expect the lensing 
signal to predominantly improve the constraints when 
combined with observations of the primary CMB signal. 
From Fig. [S] we expect largest improvements on early 
dark energy, but less so on c^ is and Cg, except when they 
both take on low values. We verify this numerically in 
Sec.lTV] 



B. CMB Lensing 



C. CMB Experiments and Forecasting 



Gravitational lensing of the CMB can improve sig- 
nificantly the CMB constraints on several cosmological 



To evaluate the future constraints on EDE models we 
consider the Planck (2|| and CMBPol [27| experiments 
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using three frequency channels for each with the experi- 
mental specifications as listed in Table U below. 



Experiment Channel[GHz] 


FWHM 


Ot [fJ-K] 


a P [fiK 


Planck 


143 


7.1' 


6.0 


11.4 


fsky = 0.85 


100 


10.0' 


6.8 


10.9 




70 


14.0' 


12.8 


18.3 


CMBPol 


150 


5.6' 


0.177 


0.250 


fsky — 0.85 


100 


8.4' 


0.151 


0.214 




70 


12.0' 


0.148 


0.209 



TABLE I: Planck and CMBPol experimental specifications. 



We consider for each frequency channel a detector noise 
of (9a) 2 where 8 is the FWHM of the beam assuming a 
Gaussian profile and a is the sensitivity. We therefore 
add to each Cg, fiducial spectrum a noise spectrum given 
by: 



(9a x ) 



2J(l+l)/ll 



(7) 



where lb = V° In 2/9 and the label X refers to either 
temperature or polarization, X = T, P. 

When CMB lensing information is also included we 
add to our dataset the lensing deflection angle power 
spectrum (and the corresponding noise spectrum). At 
sufficiently large angular scales (I < 1000), contributions 
to the deflection field will come mainly from the linear 
regime and, in harmonic space, the power spectrum of 
the deflection field reads : 



(at4 m ,) = (C? d + N{ ld )8 w S T1 



(8) 



where af m can be considered as an approximately Gaus- 



sian variable |34| . The noise power spectrum N dd reflects 



the errors in the deflection map reconstruction. We es- 
timate the lensing contribution with the quadratic es- 
timator method of Hu & Okamoto 



34] 



based on the 

correlations between five possible pairs of maps: TT, 
EE, TE, TB, EB (since the B-mode signal is domi- 
nated by lensing on small scales, the estimator BB can- 
not be used in this method) . N dd corresponds to the min- 
imal noise spectrum achievable by optimally combining 
the five quadratic estimators. Finally, the non- vanishing 
correlations between the temperature and the deflection 
maps are : 



a lm a l'm') - °; °U' d : 



Tdi 



(9) 



Following the description in [30[ we generate Cf d , Cj d 
and N dd power spectra and include these datasets in the 
analysis, both for Planck and CMBPol. 

To get a general sense of the parameter constraints and 
degeneracies, we first perform a Fisher matrix analysis. 
The Fisher matrix is defined as : 



F 



d 2 ln£ 

d-pid-p 



3 ' Po 



(10) 



where £(data|p) is the likelihood function of a set of pa- 
rameters p given some data; the partial derivatives and 



the averaging arc evaluated using the fiducial values po 
of the parameters. The Cramer-Rao inequality implies 
that (F~ 1 )a is the smallest variance in the parameter pi, 
so we can generally think of F^ 1 as the best possible co- 
variance matrix for estimates of the vector p. The one 
sigma error for each parameter is then defined as : 



°p, = y/iF- 1 )* 



(11) 



The Fisher matrix for a CMB experiment is given by 
(see [H) : 



1=2 - » ° Pl 



>.l 



*/3 



dp i ' 



(12) 



where a and (3 are running indexes over the angular 
power spectra C/. For example we include temperature 
TT, temperature-polarization TE, E mode polarization 
EE, or TT, Td, dd in the case with CMB lensing. Cov z 
is the spectra covariance matrix. We use information in 
the power spectra out to Z ma x = 3000. 



IV. RESULTS 
A. Constraints from Planck and CMBPol 

We consider a set of 9 cosmological parameters with 
the following fiducial values: the physical baryonic and 
cold dark matter densities relative to critical il^h 2 — 
0.02258 and Q c h 2 = 0.1109, the optical depth to 
reionization r = 0.088, the Hubble parameter Hq = 
71km/s/Mpc, the current dark energy equation of state 
wq = —0.90, the early dark energy density relative to 
critical fl c = 0.03, the spectral index n s = 0.963, and 
finally the effective and viscous sound speeds c 2 and c^ is . 
In order to check the stability of the result under the 
assumption of the fiducial values for c 2 and c 2 - ls we in- 
vestigate several different pairs of values. CMB lensing 
is always included except for the comparison in Tabic HU 

Using the method described above we forecast the 
constraints on wo and £l e . We find that both Planck 
and CMBPol can constrain with high accuracy those 
parameters. Planck will obtain a^ 3,11011 = 0.10 while 
CMBPol can improve this by an order of magnitude 
to cr™ BPo1 = 0.01. The density in EDE will also be 
well constrained by Planck, with a P lanck — 0.004, while 
CMBPol can improve by a factor four to ct^ MBPo1 = 0.001 
(see also Table llllj) . We find no significant dependence 
of these constraints on the choice of the fiducial values of 
the EDE perturbation parameters c 2 and c 2 is . Figure [7] 
shows the 2-dimcnsional likelihood plots in the wo-!2 e 
plane for both Planck and CMBPol experimental config- 
urations. These results are for the case c^ is = c 2 = 0.33, 
but again, there is no practically change in the contours 
for different choices of c 2 or c 2 is . 

The expected 1-a constraints on EDE perturbation pa- 
rameters c 2 and c 2 is arc presented in Table HT1 for Planck 
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FIG. 7: 68% and 95% c.l. likelihood contours for Planck (solid 
line) and CMBPol (dashed line). The ' + ' symbol represents 
the fiducial values. 



and for CMBPol experiments. We show the constraints 
obtained both with and without CMB lcnsing data. 

From the results listed in Table [TT| we can derive the 
following conclusions about estimating c 2 is and c 2 : 

• Including CMB lcnsing improves the constraints by 
- 10 - 20% (as compared to 50 - 60% for Q e and 
5 - 10% for w ). 

• CMBPol provides constraints that are generally a 
factor ~ 2 better than Planck. 

• The constraints on c 2 (or c 2 is ) depend strongly 
on the assumed value of c^ is (respectively cf), the 
general trend being that the uncertainties grow 
with the fiducial values. For example, assuming 
c 2 = 0.1, the 1-cr error on this parameter will in- 
crease by a factor ~ 5 if the fiducial model moves 
from c 2 is = 0.01 to c^ is = 1. At the same time, as- 
suming c 2 is = 0.1, the 1-cr error on this parameter 
will increase by a factor ~ 2 if the fiducial model 
moves from c 2 = 0.01 to c 2 = 1. 

• The strong correlation between c 2 and c 2 ia makes 
it difficult to precisely measure these parameters 
individually with either Planck or CMBPol (and of 
course the situation worsens as O e decreases or wq 
approaches —1). 

The correlation between EDE perturbation parameters 
can be clearly seen in Figure [51 where we plot the 68% 
and 95% c.l. 2-D likelihood contour plots in the c 2 -c 2 is 
plane. The solid lines are the constraints derived from 
Planck while the dashed lines are from CMBPol. A rea- 
sonable way of quantifying how well the sound speed and 
viscosity sound speed can be constrained is by asking at 
what significance level a non-standard value of c 2 or c 2 is 
can be distinguished from the standard (quintessence) 
value, i.e. from c 2 = 1 or c 2 is = 0. By this metric, 
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FIG. 8: 68% and 95% c.l. likelihood contours for Planck (solid 
line) and CMBPol (dashed line). In the upper panel the fidu- 
cial values are cj is = 0.1 and = 1, in the middle one c^ is = 1 
and Cg = 0.1, in the lower panel an intermediate case with 
c vis — c s — 0.33 is reported. The '+' symbol represents the 
fiducial values. Note the different scales. 



whether or not the Planck and CMBpol experiments pro- 
vide much insight of course depends on the fiducial values 
of c 2 and (? vis . For example, for the c 2 = 0.1 and <? vis = 1 
fiducial model (middle panel), Planck could rule out a 
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Fiducial Fiducial 


Planck CMBPol Planck CMBPol 


Planck CMBPol Planck CMBPol 


<£■ 

Wis 




IT 9 

c Zis 


rr 9 


(7 r 2 




IT 9 

c Zia 


rr 9 


rr ■ < 


(J 2 

C s 


0.01 


0.1 


0.019 


0.008 


0.027 


0.013 


0.016 


0.007 


0.023 


0.010 


0.1 


0.1 


0.075 


0.037 


0.093 


0.043 


0.067 


0.038 


0.082 


0.031 


0.33 


0.1 


0.17 


0.081 


0.11 


0.064 


0.16 


0.092 


0.10 


0.051 


1 


0.1 


0.52 


0.27 


0.12 


0.074 


0.42 


0.20 


0.11 


0.057 


0.33 


0.33 


0.24 


0.14 


0.16 


0.11 


0.21 


0.12 


0.15 


0.10 


0.1 


0.01 


0.094 


0.048 


0.029 


0.014 


0.084 


0.032 


0.022 


0.012 


0.1 


0.1 


0.075 


0.037 


0.093 


0.043 


0.067 


0.038 


0.082 


0.031 


0.1 


0.33 


0.098 


0.061 


0.10 


0.074 


0.092 


0.058 


0.11 


0.072 


0.1 


1 


0.19 


0.10 


0.71 


0.35 


0.17 


0.091 


0.68 


0.33 



TABLE II: Fisher analysis results at 68% c.l. for several different values of and c^, is , for Planck and for CMBPol datasets, 
with and without CMB lensing included in the analysis. 



c lis = °) at about 2ct and CMBpol 
could do this at more than 4o\ However, for the <?„ = 1 



perfect fluid (i.e. c 



and 



0.1 fiducial model (top panel), neither ex- 



periment can rule out a perfect fluid. Similarly, for the 
fiducial in the middle panel, both experiments can rule 
out c 2 s = 1 (quintessence) at very high significance, but 
not for a fiducial value of c 2 s significantly closer to unity. 

The fact that the uncertainties and ellipse shapes de- 
pend strongly on the fiducial parameter values means 
that the Fisher matrix evaluated at the fiducial model 
is not a good predictor of the shape of the likelihood 
function away from the fiducial model (and that the like- 
lihood function is thus far from Gaussian). Hence, away 
from the fiducial, the true constant likelihood contours 
could be quite different from the ones calculated using 
the Fisher matrix. This means one has to be cautious 
when making estimates as in the previous paragraph. 
For example, from the middle panel of Figure [8j we esti- 
mated that Planck would rule out c 2 is = at about 2a, 
i.e. at 95% confidence level. However, since the uncer- 
tainty in c 2 is decreases strongly as the fiducial value is 
lowered, the true significance may in this case be higher 
than 95%. However, this subtlety does not affect the 
main point made in the previous paragraph, namely that 
for a range of reasonable values of c 2 s (c 2 is ) , both Planck 
and CMBpol will be able to rule out the canonical value, 
although CMBpol with much more significance. In sec- 
tion [iVEl we check our Fisher results using an MCMC 
analysis of the true non-Gaussian likelihood and we find 
that our Fisher estimates of uncertainties and error el- 
lipses calculated are quite accurate. 



B. Including Supernovae 

Since the early dark energy component changes the 
Hubble parameter and luminosity distances, Type la su- 



pernovae (SN) information can be very useful to break 
geometrical degeneracies. 

Each SN magnitude measurement can be expressed as: 

m, = 5log w [H dL(zi,w ,Q, m ,Cl e )] + M + e t (13) 

where is the luminosity distance, M is a combina- 
tion of the SN absolute magnitude and Hubble constant, 
and e is a zero mean random term including all system- 
atic and measurement errors. Given N SN at redshifts 
Z\...Zn, we can describe the measured data as an N- 
dimensional vector m. Assumin. 
Fisher matrix is given by (see |36|) 



Gaussian errors e$, the 



F 



SN 



I, 



=- 2 ^[c 



dpi dpj 



dpi dpj ' 



(14) 



where fi = (m) is the vector of mean magnitudes and 
C = (mm T ) — fi/j, T is the covariance matrix of magni- 
tudes. The parameter vector p for the SN Fisher matrix 
includes £l m , Q c , wq, and the nuisance parameter M. 

For future SN data we consider 1800 SN out to z = 1.5 
(roughly with a cut SNAP distribution j37|) plus 300 
local (z = 0.05) SN, with an intrinsic dispersion of 0.1 
mag and a systematic error of 0.02(1 + z)/2.7 per 0.1 bin 
in z added in quadrature. New EDE parameters errors, 
reported in Table IIII1 arc estimated considering a total 
Fisher matrix : 



yrOT 



CMB 



(15) 



We see that the main improvement of adding SNe is on 
wq, reducing the Planck uncertainty by a factor of 5, and 
the CMBPol one by a factor of 2. The SN measurements 
do not reach to high enough redshift to have a good han- 
dle on Q c (the distance out to z = 2 in a model with no 
early dark energy but w a = — (l/2)dw/dlnaL = i = 5f2 
agrees nearly exactly with an EDE model [25[, and w a 
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Parameter 


Planck 


CMBPol 


uncertainty 


alone +SN 


alone +SN 




0.10 0.02 


0.010 0.005 


crn e 


0.004 0.003 


0.001 0.001 


a c 2 


0.15 0.15 


0.10 0.09 




0.21 0.20 


0.12 0.11 









TABLE III: 68% c.l. uncertainties on EDE parameters from 
Planck or CMBPol with and without SN distance information. 
The fiducial values ct — c„ is = 0.33 are used. 



Model 


Planck 


CMBPol 


fi c = 
— 


0.09 

0.004 

0.007 0.20 


0.02 

0.001 

0.003 0.07 



TABLE IV: 68% c.l. uncertainties on EDE density and neu- 
trino density from Planck and CMBPol, including marginal- 
ization over the perturbation parameters and cj is . 



cannot be determined so precisely). We also see that the 
perturbation parameters appear to be mostly uncorrc- 
latcd with any parameters to which SN distances are sen- 
sitive (indeed, they will be correlated mostly with each 
other). It is not clear what probes are best for further 
constraining c 2 and c^ is , since CMB lensing (especially at 
the level of CMBPol) already includes matter power spec- 
trum information. Perhaps three-dimensional weak lens- 
ing and galaxy statistics, or nonlinear structure, would 
supply more leverage. We leave this for future work. 



C. Including Massive Neutrinos 

In addition to considering situations where the per- 
turbation parameter constraints improve, we should also 
explore other parameters that might be degenerate with 
them, and so both weaken the constraints and be affected 
themselves by the presence of cold or stressed dark en- 
ergy. 

In particular it is interesting to study whether EDE 
could have any implication for the bounds on the neu- 
trino mass from CMB experiments. Planck and CMBPol 
are indeed expected to provide new and very stringent 
bounds on the sum of neutrino masses ^ , extremely 
competitive with respect to bounds coming from labora- 
tory experiments as KATRIN (3§ |. 

We performed a new Fisher matrix analysis adding 
to our 9-dimcnsional set of cosmological parameters the 
neutrino energy density, fi^/i 2 , with a fiducial value of 
0.001 (corresponding to ^2m u sa 0.09 eV; we quote all 
results in terms of J2 m f = 94f2 ; ,ft 2 eV). In Figure [3] we 
report the constraints from Planck and CMBPol and as 
we can see there is an anticorrelation between f2 e and 
Y^, fn v for both Planck and CMBPol experiments. This 
means that future CMB bounds on the neutrino mass can 
be affected by the presence of an EDE component (also 
see [39| ) . Numerical results are reported in Table IIVI In 
particular, we studied the impact of one component on 
the other. As we can see from the Table the presence of 
early dark energy and massive neutrinos almost doubles 
the uncertainty on both of these parameters. 

Moreover a wrong assumption of the Q c fiducial value 
(e.g. ignoring early dark energy) can bias the estimation 




0.01 0.02 0.03 0.04 0.05 

FIG. 9: 68% and 95% c.l. likelihood contours for Planck (solid 
line) and CMBPol (dashed line). The degeneracy between fi e 
and Yl m v means that the constraint on Y m v is affected by 
the inclusion of an EDE component, by a factor 2-3 as seen 
in Table El 

of other parameters and in particular of neutrino mass, 
as we now discuss. 



D. Bias from Neglecting Perturbations 

With the Fisher matrix formalism we can also evaluate 
the bias generated in parameter estimation when analyz- 
ing the datasets assuming a wrong fiducial model, e.g. 
fixing c^ is to the wrong value. 

For a Gaussian likelihood function, the bias in the i-th 
cosmological parameter, S9i, caused by the discrepancy 
between the assumed value of a parameter ipj and its true 
value, Stpj, is given by [4"oT - l42l | : 

m = -[F 66 ]uFlf^ (16) 

where F is the Fisher matrix in the space of Bi param- 
eters, and F e ^ is a Fisher submatrix with derivatives 
with respect to the assumed bias parameters tpj and the 
measured parameters 

In our case we want to study the effect of fixing c^ is = 
when an input ("true") model has c^ is = 0.33. Figure ITD1 



shows the shift obtained on the early dark energy pa- 
rameters Cg , wo, and Q e . We plot 2-dimensional contours 
showing the degeneracies at 68% and 95% confidence lev- 
els for Planck in the left panels and CMBPol in the right 
panels. The solid lines are the results obtained including 
Cy is in the parameter marginalization, while the dashed 
lines are the contours obtained when c^ is is (incorrectly) 
fixed to 0. 

As expected the constraint on cf can be affected by a 
wrong assumption on c^ is . Assuming a value of c^ is lower 
than the truth is like assuming more perturbations, so 
Cg must be biased high to compensate and reduce the 
perturbations. The resulting best fit value is ~ 1-tr away 
from the fiducial value for Planck, and ~ 2-er away for 
CMBPol. The other parameters are only mildly biased. 

When massive neutrinos are considered, f2 will play 
the major role and will strongly affect ^2 m v . In particu- 
lar we study the effect of neglecting early dark energy (i.e. 
fixing Q c = 0) when an input true model with Q c = 0.03 
is used. This assumption will shift ^2 m v from its true 
value of 0.09 eV to 0.59 eV and 0.65 eV for Planck and 
CMBPol respectively - excluding the true value by 28er 
in the latter case! 



E. Comparisons with MCMC 

Because the Fisher matrix forecasts are sometimes bi- 
ased, especially in case where there is a strong degener- 
acy between parameters, we check our previous results 
with the analysis that maps out the full likelihood func- 
tion in the cosmological parameters. The analysis uses 
the publicly available MCMC package cosmomc [43[ with 
a convergence diagnostic done through the Gelman and 
Rubin statistic. We sample the following 11-dimcnsional 
set of cosmological parameters, adopting flat priors on 
them: the baryon and cold dark matter densities a>b and 
w c , the Hubble constant Hq, the scalar spectral index ns, 
the overall normalization of the spectrum A at k = 0.05 
Mpc -1 , the optical depth to reionization, t, the current 
equation of state parameter wq, the early dark energy 
density f2 e , the dark energy sound speed logCg, the vis- 
cosity sound speed cj is , and the neutrino masses m v ■ 
We consider purely adiabatic initial conditions and we 
impose spatial flatness. We moreover only consider u>o 
values greater than —1. The fiducial model for generating 
the mock data uses the WMAP seven year best fit cos- 
mological parameters values, plus wq = —0.9, Q e = 0.03, 
c^ is = 0.33 and = 1. 

The results obtained are in good agreement with Fisher 
constraints, recovering the fiducial value at the ltr level 
for all the parameters. Moreover, the MCMC errors are 
in good agreement with the Fisher matrix error esti- 
mates, as reported in Table fVl 
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Parameter 


Fisher 


MCMC 


w 


0.10 


0.10 




0.004 


0.007 


d 


0.73 


0.74 


c 2 . 

Wis 


0.26 


0.27 



TABLE V: l-cr errors from Fisher matrix and MCMC analysis 
on EDE parameters from Planck dataset. 



V. DISCUSSION 

In this paper we have investigated future constraints 
on EDE models achievable by Planck and CMBPol ex- 
periments. We included CMB lensing as a probe, and the 
possibilities of a sound speed less than the speed of light 
and of anisotropic stresses in the clustering of the dark 
energy component parameterized with a viscosity param- 
eter c^ is . Overall, the model can be viewed as "early, cold, 
or stressed dark energy" . 

We have found that c^ is can be strongly correlated with 
the sound speed parameter cf. For this reason it will 
be difficult for these future experiments to derive signif- 
icant constraints on these sound speed parameters indi- 
vidually, although finding a deviation from the standard 
quintessence with cf = 1, c^ is = will be possible. 

We have also shown that neglecting the possibility of 
anisotropic stresses in EDE could significantly bias the 
constraints on EDE parameters. 

The results, obtained through a Fisher Matrix formal- 
ism, have been checked by a Monte Carlo Markov Chain 
analysis on Planck synthetic data. We have considered 
SN information to break geometrical degeneracies and 
we have found this significantly improves the equation of 
state parameter estimation. Finally we have investigated 
the impact of EDE on the determination of the neutrino 
mass from CMB experiments and we found it to be sig- 
nificant. In particular, neglect or misestimation of early 
dark energy density can severely bias neutrino mass con- 
straints for both Planck and CMBPol. Investigation of 
early, cold, or stressed dark energy is important not only 
to uncover further windows on the nature of dark energy 
and high energy physics, but to ensure that conclusions 
on other cosmological parameters are robust. 
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